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ABSTRACT 

The daily averages of total column amount of 
ozone taken in the period 1964-1988 at a network of 24 
Dobson stations have been analyzed. Year-round data 
as well as summerdata (May-Aug.) and winter data 
(Dec. -March) have been examined in the following 
regions: latitude bands (30°N-39°N, 40°N-52 0 N, 30°N- 
6 0°N), North America, Europe, Japan. 

Tofind year-to-year changes in the shape of the 
annual statistical distribution of total ozone (ASDTO) 
for these regions, we analyze trends in the following 
statistic characteristics of ASDTO: mean, standard 
deviation, median, 10 and 90 percentiles. Time series 
of the statistical c h a r a c t e r i s t i c s f o r the selected 
regions have been combined by averaging the 
individual stations values of these characteristics. 

The trends have been calculated by the multiple 
regression model adjusted for: the 11-year solar cycle, 
the Southern Oscillations effects, and for serial 
correlations. We have found that: 

a) in all regions (excluding Japan, North America), 
the shape of ASDTO has been drifting towards low 
ozone values. The drift seems to be not accompanied 
with a transformation in Ihe shape of ASDTO. 

The drift speed (therate of decrease in the annual 
means oftotal ozone) is of order 1-3 % per decade (in 
the period 1970-1988). 

b) the drift speed depends on the regions, the lowest 
one is in 3 0°N-3 9°N band. 

b) in Japan, the interannual changes in the shape of 
ASDTO have not been revealed. 

c) in North America, the drift of the year-round 
ASDTO (the year-round ASDTO comprises all the daily 
means oftotal ozone in a given year) has been 
accompanied with a transformation in the shape. The 
shape of the year round ASDTO becomes narrower. 

d) in all regions, except Japan and the band 30P-3 9°N, 
the winter ASDTO (the winter ASDTO comprises the 
date taken in the period December in a given year 
through March next year) moves faster towards low 
ozone values than the summer ASDTO(the summer 
ASDTO comprises the data taken in the period May- 
through August in a given year). 

1. INTRODUCTION 

Ozone trend detection has been a subject of 
interest during the past decade (e.g. Reinsel et a l ., 
1981; Angel 1 et a/., 1983; Hiller al., 1986; Oehlert, 1986; ’ 


IOTP, 1990). Implications for UV radiation changes at 
the earth's surface have stimulated muchof the 
concern about ozone downward trend. 

Almost all studies of the long-term variations in the 
atmospheric ozone were carried using the monthly 
(or yearly) averagesof ozone. In this paper, we 
analyze interannual changes in valuesof a few 
statistic characteristics of the annual statistical 
distribution of total ozone (i.e. mean value, standard 
deviation, median, 10 and 90 percentiles) to disclose 
the long-term variations in the shape of ASDTO in 
selected Northern Hemisphere regions 

2. DATA 

Dobson total ozone data have been analyzed 
from a network of 24 stations for the period 1964- 
1988. The names of the stations, locations and the 
data periods are listed below: 


STATION 

COUNTRY 

LAT. 

LON. 

PERIOD 

1 Lerwick, 

U.K. 

60°N, 

1°W ; 

1/64-1 1/88 

2 Leningrad, 

Russia 

60°N, 

30°E; 

8/68-12/88 

3 Churchill, 

Canada 

59°N, 

94° W; 

1/65-12/88 

4 Edmonton, 

Canada 

54°N, 

1 1 4°W; 

1/64-12/88 

5 Belsk, 

Poland 

52°N, 

21°E; 

1/64-12/88 

6 Bracknell, 

U.K. 

5 1°N, 

1°W; 

1/69-12/88 

7 Ucclc, 

Belgium 

5 1°N, 

4°E; 

2/71-12/88 

8 Hradec K., 

Czechoslov. 

50°N, 

16°E; 

1, 64-12/88 

9 Hohenpciss 

..Germany- 

48°N, 

1 1°E: 

1/67-12/88 

10 Caribou, 

U.S.A. 

47°N, 

68°W ; 

1/64- 1 2/88 

1 1 Arosa, 

Switzerland 

47°N, 

10°E; 

1/64-12/88 

12 Bismarck, 

U.S.A. 

47°N, 

1 0 1° W; 

1,64-12/88 

13 Toronto, 

Canada 

44°N, 

79°W : 

1/64-12 88 

14 Sapporo, 

Japan 

4 3°N , 

1 41°E; 

1/64-12,88 

15 Rome, 

Italy 

42°N, 

12°E; 

1/64-12,88 

16 Boulder, 

U.S.A. 

40°N, 

1 05°W : 

1. 64-12 88 

17 Cagliari, 

Italy 

39°N, 

9°E ; 

1,64-12,88 

18 Wallops Is., 

U.S.A. 

38°N, 

76°W ; 

1, 70-12. 88 

19 Nashville, 

U.S.A. 

36°N, 

8 7° W ; 

1 , 64-12,88 

20 Tateno, 

Japan 

3 6° N , 

14(^E; 

1,64-12,88 

21 Srinagar, 

India 

34°N, 

75°E; 

2, 64-12. 8 8 

22 Kagoshima, 

Japan 

32°N, 

1 3 1 ” E : 

1 64- 1 2 88 

23 Ouctta, 

Pakistan 

30°N, 

6 7’E; 

8, 69- 12. 8 8 

24 Cairo, 

Egypt 

30°N, 

3 1°E: 

1 1 74-12,88 
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In this paper, we analyze the daily mean total ozone 
record published in "Ozone Data for the World" (ODW) 
journals. However, International Trend Panel-1988 
(IOTP, 1 990] found strong incongruities in the 
published data of many ground-based stations due to 
instrument calibrations. 

For 31 Dobson stations, a provisionally revised data set 
(195 7-1986) was prepared by the Panel, with the 
corrections being applied to the monthly averages of 
total ozone. 

We examine results of theozone observations taken at 
the stations selected by the Panel. 

Making use of the Panel's re-evaluated total ozone 
monthly means we convert daily averages of total 
ozone (calculated before January 1987) to the 
provisionally corrected ones in the following way: 

. , f Ml (j,k,l) 

°3 n e w (* J* = 03 old ( 1 >J > *’ *) M2(j,k y l) 

w h e r e : Oj n e w ( 1 ,j , k , /), Oy 0 / </ ( I , j , k , l) - d a i 1 y 
mean of total ozone o n i -t h d a y , j - 1 h month, k-th 
year and at/ -th Do b so n station, provisionally 
revised and uncorrected, respectively ,Ml(j,k,l)- 
rc vised monthly mean of total ozone [IOTP, 1990], 
M2(j,k,l)- monthly mean of totalozone derived 
from ODW data. 

Provisionally revised daily averages of total ozone 
have been transformed to the departures from the 
long-term monthly means ( 1 964- 1 988) expressed in 
percents of the long-term monthly standard 
deviation. 

Using Oj „ e w (i,j,k, I) and 0 3 0 i j (i,j,k,l) (since 
January 198 7) we calculate a few characteristics 
describing the shape of ASDTO at station / in year*. 
The following ones have been considered: 
mean value(Af£L4;V), standard deviation (S D) , median 
{MED), 90 percentile ( P90 ), 10 percentile (P 1 0) P 1 0 
and P 9 0 estimate locations ofthe low values tail and 
the high values tail of ASDTO, respectively. 

For selected northern hemisphere regions, the 
statistic characteristics of ASDTO in a given year have 
been estimated averaging the individual stations 
values of the statistic characteristics in that year. 

The following regions have been considered : latitude 
bands ( 3 0°N - 3 9°N , 40°N-52°N, 30°N-60°N), North 
America, Europe, Japan. 

The statistic characteristics of ASDTO have been 
calculated from all daily averages of total ozonevalucs 
in a given year. To find seasonal differences inthe 
long-term changes of ASDTO, summer (May through 
Aug.) and winter (Dec. through March) values ofthe 
statistic characteristics have been examined also. 

3. TREND MODEL 

Trends in the statistic characteristics of ASDTO 
have been found by m u 1 1 i p 1 e reg re s si o n model. The 
model has been adjusted for: the 11-year solar activity 

cycle, the long-term stratospheric circulation 

fluctuations (Southern Oscillations) effects, andthe 
serial correlations. 

Variations in total ozone related to the QBO have not 
been parameterized because the r e 1 a t i v e 1 y h i g h 
frequency of the QBO (there are about 12cycles in the 
period 1964-1988) makes its inclusion in the ozone 
trend model (or exclusion) have a negligibleeffect on 
the trend estimates (IOTP, 1990). 


Trends in the statistic characteristics of ASDTO 
havebeen modelled by means of the ramp trend 
model, i.e. we assume that the trend in ozone has 
formed after 1969: 

y (t) = a+0Sun(t)+rSoi(t)+N(t) y t ^1969 

y(t) = aj + Tt + pSun(t) + ySoi(t)+ N(t), t >1969 

where: a, aj f fi* Y ~ constants, T - rate of linear 
changes (trend), S u n ( t )- annual mean ofZurich 
sunspot number in year/, S o i ( t )- Southern 
Oscillation index (normalized press u re difference 
between Tahiti and Darwin). The noise term N ( t ) i s 
modeled as a first order autoregressive process, i.e. 
N (t) = 6N ( t - 1 ) +e ( t ), where is a noise with 

normal distribution and zero mean. 

In our model, we assume that declining tendency in 
total ozone has begun around 1970. This hypothesis 
was widely used in the trend modeling in thepast 
decade (e.g. Reinsel et al. t 1981*, IOTP, 1990). 

We run the trendmodeJ using as jr f f J the composite 
time series of: MEAN, S D, MED , P 9 0, and P/0, 
combined from individual station values of these 
characteristics. Constants a, ctl,fi, Y> ^ T are derived 
by the least square fit. 

For example, the comparison between observ ed and 
modeled lime ser ies o f MEAN f 0 r 10 European 
Dobson station (see section 2 for the names of these 
stations) is presented in Fig. l.The correspondence 
between the long-term variations in the observed 
and modeled time series allows to state that the model 
reproduces adequately the long-term variations 
in total ozone. 



Fig.l Observed and modelled departures of annual 
total ozone means from its long-term (1964-1988) 
mean (expressed inpercent of the long-term mean) 
in the period 1 9 64 - 1 9 8 8 for European Dobson 
stations (see section 2 for the names of stations) 

4. RESULTS 

In Fig. 2 we present an example of analyzed 
time series set. The annual values of the statistic 
characteristics h a v e been combined from ten 
European Dobson stations. It is seen that the pattern 
o f SD time series is almost c o n s l a n t d u r i n g 
observation’s period ( 1 964- 1 9 88), b u t other time 
series {ME A N, P 9 0, P 1 0 ) show similar decreasing 
tendency. Moreover, variations in time series of 
MEAN, P90 and PI 0 seem to be highly correlated. 
The correlation coefficients be t we e n : ME A N and 
P9 0, MEAN andP/0, P^andP/Oare equal 0.93, 
0.93, 0.75, respectively. 
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These high correlations, similar negative trends in 
MEAN f PI0,P90 , and no trend in 5 Diet us suggest 
that year-to-year variations in the shape ofASDTO are 
caused by a drift ofthe shape (towards low ozone 
values) rather than by a transformation of the shape 
(for example a change of the distribution type). 

For the European stations, we have found that ASDTO 
is well approximated by the log-normal distribution. 
The long-term changes in the shape of ASDTO canbe 
mathematically de scribed as a translation of the 
distribution towards low ozone values because almost 
similar negative trendshave been found in time 
series of MEAN, P 1 0 and P 9 0, while no trend in SD. 

For allanalyzed northern hemisphere regions, 
the high correlations between A/K/tiV, P/Oandi’PO 
have been detected. Then, the long-term changes in 
the shape of ASDTO have been deduced from the trend 
analyses of the following c o m p o s i t e t i m e series: 
MEAN , SD, MEDIAN , P 1 0, P90. The results are 
presented in Tab. 1 . 



Fig. 2. Time series (1964-1988) of statistic 
characteristics (annual me an- MEAN, standard 
deviation- SD, 90 percentile- P 9 0 ,10 percentile -P 1 0) 
of the distribution of total ozone combined from the 
individual station values of these characteristics. The 
analyzed time series have been obtained from the data 
taken at 10 European Dobson stations. Values of the 
characteristics are expressed in standard deviation 
unit. Dashed lines represent linear regression lines. 

The main finding fromTab.l is a drift of the 
shape ofASDTO towards low ozone values (only the 
results for Japan showthat the shape remains 
constant in the period 1 970- 1 988). We have found that 
theannual rates of the changes in selected statistical 
characteristics (MEAN, MEDIAN, P 9 0, P 1 0 ) ofASDTO 
are almost the same (within 2o band) and the trends 
in SD are not statistically significant (except North 
America where the trend is negative). Therefore, it 
seems that the long-term changes inthe shape of 
ASDTO in all the analyzed regions (except Japan, 
North America) are expressed simply as a translation 
of the shape towards low ozone values. The speed of 
this translation (the rate of decrease in the annual 
means of total ozone) depends on latitude (the lowest 
is found in the latitude band 3 0°N-3 9°N) The speed 
varies with season also. For the winter ASDTO, which 
is built from the all daily means of total ozone in the 
period December- M arch, the speed is greaterthan the 
speed for the summer ASDTO (Maythrough August 
daily means of total ozone have been usedto prepare 
the summer ASDTO). 


I n North America, for the year round data, 
stat isti cal ly s ig n i f i c an t trend in SD time series has 
been revealed. Therefore, the drift of the shape of 
the year-round ASDTO (the year-round ASDTO 
comprises all the daily means of total ozonein a 
given year) seems to be accompanied with a 
transformation of the shape. The shapeof the year 
round ASDTO becomes narrower. There is a 
likelihood that thenegative tendency in the annual 
s t an d a r d de v i at i on s over North America, calculated 
from all the daily total ozone means in a given year, 
is forced by the large decline in the high value tail 
of ASDTO (see Table 1, section for North America, 
where a large difference between the trends in P 9 0 
and P10 can be seen). 

Tab. 1 . Long-term trends (year-round, winter a nd 
summer trends in the period 1970-1988) inthe 
characteristics of the annual statistical distribution 
of total ozone (in % per decade) for selected 
northern hemisphere regions: latitude bands (3(PN- 
3 9°N , 4 0°N-52°N , 3 0°N- 60°N ) , North America , 

Europe, Japan. 

Forthe period 1964-1986 provisionally revised daily 
averages of total ozone and for the period 1987-1988 
the data published in Ozone Data for the World 
journals have been used. Errors (on 2o level) are 
shown in parentheses. Statistical significant results 
are underlined 


EUROPE 

Trends_ 



_Year-round_ 

_Winter_ 

_S u m m e r _ 

MEAN 

.80) 

-2.50 (2.89) 

-1.11 (1.83) 

SD 

-1.64 (5.66) 

2.54 (7.13) 

-3.10 (7.20) 

MED 

-1.98 (1.82) 

-1.96 (2.99) 

-1.20 (1.79) 

P90 

-2.45 (2.32) 

-2.17 (3.65) 

-1.63 (2.29) 

P10 

-1.99(1 .60) 

-3.33 (2.37) 

-0.70 (1.72) 


NORTH AME 

RICA 



Trends 



Year-round 

_Winter_ 

_Summer_ 

MEAN 

-2.60(1.62) 

-3.83 (2.14) 

,2.03,(1.4 7) 

SD 

-4.91 (3.38) 

-1.44 (6.12) 

-5.29 (5.87) 

MED 

-2.70 (1.59) 

-UL5_(2.12) 

-2.17 (1.46) 

P90 

-3.26 (1.98) 

4.11.(2.80) 

_ -2.5.0 (1.98) 

P10 

-1.95 (1.40) 

-_2JfT(2.16) 

1.3.7 (1.37) 



_Year-round 

JAPAN 

Trends 

_Winter_ 

Summer 

MEAN 

-0.70 (1.19) 

-0.73 (2.01) 

0.05 (1.68) 

SD 

0.59 (4.30) 

-2.67 (9.14) 

0.16 (4.61) 

MED 

0.16 (1.25) 

-0.88 (2.12) 

0.33 (1.78) 

P90 

0.06 (1.29) 

-0.44 (2.80) 

0.33 (2.10) 

P10 

0.14 (1.15) 

-0.11 (1.89) 

-0.07 (1.33) 
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BAND: 40°N - 52°N 


Trends 



_Year-round_ 

_ Winter_ 

_Summer_ 

MEAN 

--1.9Q (1.35) 

-2.33 (2.06) 

-1.34 (1.23) 

SD 

*2.83 (5.12) 

-3.18 (7.43) 

-3.53 (6.24) 

MED 

-1.82 (1.34) 

-2.11 (2.11) 

-l._44_(l .20) 

P90 

-2.27(1.77) 

-2.21 (2.79) 

-1.50 (1.65) 

P10 

-_1^4_(1 12) 

-3.01 (1.71) 

-0.90 (1.23) 


BAND: 30° N - 39°N 




Trends 



Year-round 

Win ter_ 

_Summer_ 

MEAN 

-1.34(1.13) 

-1.52 (2.14) 

-1.50(1.50) 

SD 

-1.68 (4.36) 

2.00 (5.93) 

0.96 (5.62) 

MED 

-1.34 (1.19) 

-1.51 (2.23) 

,1.5241.52) 

P90 

-1.29(1.27) 

-0.84 (1.99) 

-1.51 (1.66) 

P10 

-1.44 (1.00) 

-1.45 (2.40) 

-1.46 (1.30) 


BAND: 3 0° N - 6 0° N 



Trends 



Year-round 

l_ _Winter_ 

_Summer_ 

MEAN 

-1.66(1.18) 

-2.30(1.54) 

tA^22 (1.21) 

SD 

-1.62 (3.25) 

-0.45 (4.66) 

-2.24 (3.25) 

MED 

-1.65 (1.18) 

-2.14(1.62) 

-T-.2.6 (1.24) 

P90 

-1.87(1.44) 

r2*lfi(1.84) 

-1.39 (1.44) 

P10 

-1.52 (0.96) 

-2.46 (1.34) 

-0.98 (1.01) 

5. Discussion and Conclusion 



We have found that thelong 

term change 


the shape of ASDTO, ca 1 c u 1 ate d for selected northern 
hemisphere geographical regions, have been 
expressed as a drift of the shape towards low ozone 
values Some regional and seasonal peculiarities in 
the drift specd(the rate of decrease in annual 
means of total ozone) have been revealed. 

The errors of the trend estimates in the statistical 
characteristics of ASDTO have been calculated too 
large to make conclusive statement about details of a 
transformation of ASDTO in the period 1964-1988. 
However, in all the analyzed regions (excluding 
Japan, North America), the high correlations 
between time series of MEAN, P 9 0 and P l 0 and 
almost similar trends inthese characteristics (while 
no trends in 5/)) let us suggest that the shape of 
ASDTO has been drifting towards low values without 
a transformation of theshape. Results for the North 
America stations show that the drift of the shapehas 
been accompanied with a transformation of the 
shape. The rate of ozone loss found in the high 
values tail of ASDTO(P90) is greater than ozone loss 
\ \\ M E AN, and ozone loss in the low values tail of 
PDFTO (P 1 0) is lower than the trend in MEAN. As 
the result of these changes, 50 should show 
decreasing tendency. The trend analysis ofS'Dtime 
series reveals this tendency. It means thatthe shape 
of ASDTO becomes narrowerduring its drifting 
towards low ozone values. 

The long-term trends (in theperiod 1970-1990) 
calculated from the monthly averages of Dobson 


total ozone values (Stolarski a/., 1992) are broadly 
consistent with the trends InMEAN derived from 
the interannual variations of ASDTO, i.e. there is a 
significant winter loss and a smaller summer loss i n 
total ozone over the Northern H e m i s p h e r e , t h i s 
difference is more pronounced at highand mid 
latitudes than at low latitudes. 

The results of ourtrend model (for the 
variable MEAN) and the newest trend estimates in 
monthly means ofDobson total ozone by Stolarskicf 
al. [1992] are almostthe same. The differences 
between the ozone trends are lowerthan 0.6 % per 
decade for North America, Europe, and the northern 
temperate regions. For example, Stolarski et al. 
[1992] found the following trends over 26°N-64°N 
band in the period 1970-1990: year-round *1.8% per 
decade, winter - 2.7% per decade, summer -1.3% per 
decade. 

The errors in ourtrend estimates are larger, 
because relatively small number of data points 
(annual means) hasbeen used in the trend 
detection. 

In the last decade, the rate ofozone loss 
seems to be greater than that in the period 1970- 
1990, Stolarski et al. [1992]. We have tried to prove 
this hypothesis comparing the trends in the periods 
1970-1986 and 1978-1988. The differences between 
the trends in the statistical characteristics of ASDTO 
calculated for these periods havebeen found 
statistically insignificant (limited numberof the 
data points has been used to estimate trends, 
especially inthe period 1 978- 1 988). Then, the 
analysis of the interannual variations inASDTO for 
the detection short-term trends (in the period of one 
decade)is not able to provide a support for the 
hypothesis probably due to large errors of the trend 
estimates. 
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